Adiabatic perturbation theory: from Landau-Zener problem to quenching through a 

quantum critical point 
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We discuss the application of the adiabatic perturbation theory to analyze the dynamics in various 
systems in the limit of slow parametric changes of the Hamiltonian. We first consider a two- level 
system and give an elementary derivation of the asymptotics of the transition probability when the 
tuning parameter slowly changes in the finite range. Then we apply this perturbation theory to 
many-particle systems with low energy spectrum characterized by quasiparticle excitations. Within 
this approach we derive the scaling of various quantities such as the density of generated defects, 
entropy and energy. We discuss the applications of this approach to a specific situation where the 
system crosses a quantum critical point. We also show the connection between adiabatic and sudden 
quenches near a quantum phase transitions and discuss the effects of quasiparticle statistics on slow 
and sudden quenches at finite temperatures. 



I. INTRODUCTION 



The dynamics in closed systems has recently attracted 
a lot of theoretical interest largely following the ex- 
perimental developments in cold atoms systems (see 
E.g. Ref. [l| for a review). Several spectacular ex- 
periments already explored different aspects of the non- 
equilibrium dynamics in interacting many-particle sys- 
tems @, H, 11 H, IE 0, H- Recent theoretical works 
in this context have focused on various topics, for in- 
stance: the connection of dynamics and thermodynam- 
ics [EEEEII,E3], the dynamics following a sudden quench 
in low dimensional systems [H EE EI EE EE E3, EE El, 
HE HI S3, EE E3, EH EL the adiabatic dynamics near 
quantum criticalpointsjpj^ijjl IM IM IS 11 11 IM 
M, HE H HE EE El, ElElEj Although there is so 
far very limited understanding of the generic aspects of 
the non-equilibrium quantum dynamics, it has been rec- 
ognized that such issues as integrability, dimensionality, 
universality (near critical points) can be explored to un- 
derstand the non-equilibrium behavior of many-particle 
systems in various specific situations. 

The aim of this paper is to address some generic as- 
pects of nearly adiabatic dynamics in many-particle sys- 
tems. In particular, we will discuss in details the scaling 
of the density of generated quasiparticles n GX , entropy 
Sd, and heating Q (non-adiabatic part of the energy), 
with the quenching rate. It has been already under- 
stood that this scaling is universal near quantum critical 
points [27L [28L [30l| . an d more generally in low-dimensional 
gapless systems [31j. The universality comes from the 
fact that if the system is initially prepared in the ground 
state or in a state with small temperature, then under 
slow perturbations very few transitions happen and the 
system effectively explores only the low energy part of the 
spectrum. This low energy part can be described by a 
small number of parameters characterizing some effective 
low energy theory (typically field theory). The situation 
can become different, however, in high dimensional sys- 
tems [3lj ] . This is mainly due to a typically small density 
of low energy states, which for e.g. for free quasiparti- 



cles scales as p(e) oc e d / z_1 , where d is the dimensionality 
and z is the dynamical exponent determining the scaling 
of energy e with momentum k at small k: e(k) ~ k z . 
As a result the transitions to high energy states domi- 
nate the dynamics and the universality is lost. A simi- 
lar situation happens for sudden quenches near quantum 
critical points [451 ] . In low dimensions the excitations of 
low energy quasiparticles determine the (universal) scal- 
ing of various thermodynamic quantities. However, in 
high dimensions the transitions to the high energy states 
following the quench become predominant. In this case 
one can use the ordinary perturbation theory or linear 
response which predicts that n cx ,Sd, and Q become an- 
alytic (quadratic) functions of the rate for slow quenches 
and of the quench amplitude for sudden quenches. 

In this paper we will explain in details how the tran- 
sition between quadratic and universal regimes can be 
understood as a result of breakdown of the linear re- 
sponse. More specifically we will illustrate how exactly 
the crossover between different scaling regimes occurs in 
the situations where the system can be well described 
by quasiparticle excitations. We will concentrate on the 
slow, linear in time, quenches and briefly mention the sit- 
uation with fast quenches in the end. First we will discuss 
the adiabatic perturbation theory and its implications for 
many-particle systems. Then using this theory we will 
analyze a simple driven two-level system (Landau-Zener 
problem [46l 1471]) where the coupling linearly changes in 
time in the finite range. We will show how the quadratic 
scaling of the transition probability with the rate emerges 
from this perturbation theory. Then we will consider a 
more complicated situation where the system consists of 
free gapless quasiparticle excitations. We will show that 
in low dimensions, d < 2z, the scaling of the density 
of excitations and entropy is universal n cx ,Sd ~ |<5| d ^ z , 
while in high dimensions the quadratic scaling is restored. 
The quadratic scaling can be understood as the result of 
multiple Landau-Zener transitions to high energy quasi- 
particle states. We will illustrate our argument with 
a specific model of coupled harmonic oscillators. Next 
we will consider a more complicated situation where the 
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system is quenched through a second order quantum 
phase transition. We will show how the universal scal- 
ing n cx ~ \8\ dv l( zv+1 ) [13, Hi| {v is the critical exponent 
for the correlation length) emerges from combining adi- 
abatic perturbation theory and universal scaling form of 
energies and matrix elements near the quantum critical 
point. We will also show how this scaling law breaks 
down and is substituted by a simple quadratic relation 
n cx ~ 5 2 when the exponent dv/(zv+l) exceeds two. We 
will illustrate these results using specific exactly solvable 
models. We will also discuss the connection between adi- 
abatic and sudden quenches near the quantum critical 
point. In particular, we will show that in low dimen- 
sions, dv < 2, the density of excited quasiparticles for 
a slow quench can be understood as a result of a sud- 
den quench, if one correctly identifies the quench ampli- 
tude A* (for the sudden quench) with the quench rate 
S (for the slow quench): A* ~ |^| 1 /( 2:iy + 1 ) , This analogy 
is very similar to the Kibble-Zurek picture of topologi- 
cal defect formation for quenches through classical phase 
transitions 0|4§|, where one assumes that below certain 
energy (temperature) scale topological excitations essen- 
tially freeze. However at higher dimensions, dv > 2, this 
analogy becomes misleading since for sudden quenches 
scaling of n ox is no longer determined by low energy ex- 
citations. This work mostly focuses on the situation in 
which the system is initially in the ground state. In the 
end of the paper we will discuss what happens if the 
system is initially prepared at finite temperature. We 
will argue that the statistics of low energy quasiparti- 
cles strongly affects the response of the system to fast or 
slow quenches, enhancing the non-adiabatic effects (com- 
pared to the zero-temperature case) in the bosonic case 
and suppressing them in the fermionic case. We will dis- 
cuss the corrections to the universal scaling laws if the 
low energy quasiparticles are described by either bosonic 
or fermionic statistics. 



II. ADIABATIC PERTURBATION THEORY 

We consider the following setup: the system is de- 
scribed by the Hamiltonian H(t) = Ho + X(t)V, where Tio 
is the stationary part and X(t)V is the time-dependent 
part of the Hamiltonian. Our purpose is to character- 
ize the dynamics of this system resulting from the time- 
dependent perturbation. We consider the situation in 
which the system is in a pure state. More general situ- 
ations, where the state is mixed, can be addressed simi- 
larly by either solving von Neumann's equation or aver- 
aging solutions of the Schrodinger equation with respect 
to the initial density matrix. We assume that X(t) is a 
linear function of time: 

( Xi t<0 
X(t) = l Xi + tS(Xf-Xi) 0<t<l/d . (1) 
{ X f t> 1/5 



Here 6 is the rate of change of the parameter X(t): <5 > 
corresponds to the adiabatic limit, while 5 — > oo corre- 
sponds to a sudden quench. In principle the values Xi 
and Xf can be arbitrarily far from each other, therefore 
we can not rely on the conventional perturbation theory 
in the difference between couplings |A/ — Aj|. 

In the limit of slow parametric changes, we can use 6 
as a small parameter and find an approximate solution 
of the Schrodinger equation: 

idt\1>)=H(t)W>), (2) 

where \tf>) is the wave function. Here we use the conven- 
tion that h = 1 (this can be always achieved by rescaling 
either energy or time units). Our analysis will be simi- 
lar to the one in Ref. [5fJ, nevertheless for completeness 
we will present here the details of the derivation. It is 
convenient to rewrite the Schrodinger equation ^ in the 
adiabatic (instantaneous) basis: 

M*)>=X)a„(t)K(t)>, (3) 

n 

where \(f> n (t)} are instantaneous eigenstates of the Hamil- 
tonian H(t): 

H{i)\<t> n (t)) = E n (t)\4> n (t)) (4) 

corresponding to the instantaneous eigenvalues E n {t). 
These eigenstates implicitly depend on time through 
the coupling X(t). Substituting this expansion into the 
Schrodinger equation and multiplying it by (4> m \ (to 
shorten the notations we drop the time label t in | </>„)) 
we find: 

id t a m (t) + i^a n (t)((/) m \dt\(f>n} = E m (t)a m {t). (5) 

n 

We then perform a gauge transformation: 

a n (t) = a„(t) exp [-iQ n (t)] , (6) 

where 

e„(t) = J E n (r)dT. (7) 

The lower limit of integration in the expression for ©„ (t) 
is arbitrary. We chose it to be equal to U for convenience. 
In consequence the Schrodinger equation becomes 

a n {t) = - ^a m (t){n\d t \m) exp [i{Q n (t) - 6 m (t))] . 

m 

which can also be rewritten as an integral equation: 
a n (t) = - / t dt'^a m (i')(n|a t »e J ( e "(*')- e '"( t '». 



3 



If the energy levels E n (r) and E m (r) are not degenerate, 
the matrix element (n\d t \m) can be written as 

(n\d t \m) = - jyjf^-Att)- <n|y|m> 



E n (t)-E m (t) E n {t)-E m (ty 

(10) 

where we emphasize that the eigenstates \n) and eigenen- 
ergies E n (t) are instantaneous. If X(t) is a monotonic 
function of time then in Eq. ([9]) one can change variables 
from t to X{t) and derive: 

a„(A) = - / A dA'^ am (A')(n|av|m)e l ( e "( A ')- e '»( V )) ; 

(11) 



where 



e n (A)= (\x>^± 

J \i X 



(12) 



Formally exact Eqs. (O and (TIT]) can not be solved in 
the general case. However, they allow for a systematic 
expansion of the solution in the small parameter A. In- 
deed, in the limit A — > all the transition probabilities 
are suppressed because the phase factors are strongly os- 
cillating functions of A. The only exception to this state- 
ment occurs for degenerate energy levels (5lj . which we 
do not consider in this work. In the leading order in A 
only the term with to = n should be retained in the sums 
in Eqs. (JHJ) and (TTTj) . This term results in the emergence 
of the Berry phase [52j : 

ft pMt) 
*»(*) = -il db'{n\d t ,\n) = -i dX'{n\dy\n), 

(13) 

so that 



a n (t) fa a n (0)exp[-«$„(t)]. 



(14) 



In many situations, when we deal with real Hamilto- 
nians, the Berry phase is identically equal to zero. In 
general, the Berry phase can be incorporated into our 
formalism by doing a unitary transformation a n (t) — » 
a n (t) exp[-i$ n (*)] and changing 6„ -» 9„ + $„ in 
Eqs. © and (TTT|) . 



We now compute the first order correction to the wave 
function assuming for simplicity that initially the system 
is in the pure state n = 0, so that ao(0) = 1 and a n (0) = 
for n ^ 0. In the leading order in A we can keep only 
one term with to = in the sums in Eqs. ^ and (|11|) 
and derive 



a n (t) 



or alternatively 



dt'{n\d t ,\0)e^ e ^- e ^">\ (15) 



a„(A)«- / dA>|d A '|0)e l ( e »( A V e °( A '». 



(16) 



A; 



The transition probability from the level \(f>o) to the level 
\<fi n ) as a result of the process is determined by |a n (A/)| 2 . 

The expression (|16p can be further simplified in the 
case where the initial coupling Ai is large and negative 
and the final coupling A/ is large and positive, employing 
the stationary-phase approximation. The complex roots 
of the equation E n (X*) — E (\*) = define the station- 
ary point. Consequently the dominant contribution to 
the transition probability is determined by the negative 
imaginary part of the phase difference n — 0o evaluated 
at these roots 51]: 

\a n \ 2 cx exp[-23(6 n (A*) - 6 (A*))]. (17) 

In particular, for linearly changing coupling X(t) = 5t, 
we obtain 



\a n \ oc exp 



■r 



A* 



[E n (X') - E (X')]dX' ) (18) 



and the transition probability exponentially vanishes as 
<5->0. 

However there are many cases where the coupling Xi 
or Xf or both are finite. In this case Eq. (fT7|) is no 
longer valid and the asymptotic values of the integrals 
in Eqs. (fTS"]) and (|16[) are determined by the initial and 
final times of evolution. Using the standard rules for 
evaluating the integrals of fast oscillating functions we 
find: 



a„(t f ) 



iX 



{<p n \d t \(f)o) 



d {cf>n\d t \(j)o) 



o i(e„(t)-e (t)) 



E n (t)-E (t) E n (t) - E (t) dt E n {t) - E (t) 

((f>n\d\\<Po) _^ (^nl^Al^o) ; 2 1 d (<l>n\d\Wl> 



U 



x z 



E n {X)-E (X) (E n (X) - E (X)) 2 E n (X) - E {X) dX E n (X) - E (X) 

I 



,i(e„(A)-e (A)) 



(19) 



Ai 



In the following analysis we will retain only the first non- 



vanishing term in A 



= (5. The terms proportional to 
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higher powers of the expansion parameter, such as A, 
(A) 2 , as well as non-analytic terms similar to Eq. (jTTJ) , 
will be neglected assuming sufficiently small 8 — > 0. 



The probability of the transition to the n-th level is 
approximated by: 



*n(A/)| 2 



\{K\dx t \M? 



{E n {\) - E {\)Y (E n (X f ) - E (\ f )¥ 



28 2 



(c^|d A J</>o) ((/>n\dx f \(f>o) 
E n (Xi) — Eo(Xi) E n (Xf) — Eo(Xf) 



cos [A6 n0 ] 



where A9 n0 = e il (A / )-e (A / )-6„(A l ) + e (A l ) is the 
phase difference between the states \4> n ) and |0 O ) accu- 
mulated during the time evolution. This phase difference 
is usually very large and thus the last term in Eq. ([20]) 
is a highly oscillating function, which can be typically 
dropped because of the statistical or time averaging. 



A. Application to the Landau-Zener problem. 

We apply the general formalism presented above to 
the notorious Landau-Zener (LZ) problem [4(||47|]. The 
Hamiltonian to study this problem is given by a 2 x 2 
matrix: 



H = Xa z 



(21) 



which is conveniently expressed through the Pauli matri- 
ces: 



(22) 





'1 " 




"0 f" 


Vz = 


-f 


; Ox = 


f 



This system has two eigenstates: 



_ ( sin(0/2) 
-cos(0/2) 



cos(0/2) A 
sin(0/2) ) ' 



(23) 



where tan0 = g/X, with corresponding energies E± = 

±^/F+F. 

We assume that the coupling A linearly changes in 
time: X — St. The system is initially, at t = ti, pre- 
pared in the ground state and the process continues until 
t — tf. In the limit ti — > — oo and t f — > oo the probability 
to occupy the excited state |+) [H, S3 is a non- analytic 
function of 5: 



\ a +\ 



exp 



5 



(24) 



However, in the general case, where ti oi tf are finite, 
the probability contains both non-analytic and analytic 
contributions in 5. 

In principle the LZ problem can be solved exactly for 
arbitrary t i and t / [53|, [54| (Appendix [5} , however the 
general solution is quite cumbersome. Here we illustrate 
how the asymptotical behavior of the transition probabil- 
ity at small S can be recovered employing the adiabatic 




FIG. 1: The contour of integration in Eq. (|26p in the com- 
plex f-plane. Integration over the real axis is given by several 
contributions: the one along CV around the point t = i that 
gives the first term in Eq . fl27[) . C+ and C- from the two sides 
of the branch represented by second term in Eq. (|27|l . and also 
Cr and Cr that vanish in the limit R — > oo. 

perturbation theory. The only non-zero matrix element, 
which enters Eq. (|I6p . is 



-\d t 



0/2 



I Xg 



2 A 2 



(25) 



We first apply Eq. (|T5|) to the case t 



-oo and t 



oo which corresponds to the classic LZ problem. 
Eq. (fl6|) gives: 



Then 



oo 



dt- 



Sg 



(sty 



■ exp 



2i / dr^/iSr) 2 + g< 



( 26 ) 

The asymptotic behavior of this expression at small 
8 <C g 2 can be derived by studying the analytic proper- 
ties of the integrand in the complex plane of the variable 
t = t'+it". We notice that the phase factor has a branch- 
cut singularity in the upper and lower-half planes along 
the imaginary axis, that starts at t" = ±1 and goes to 
infinity at t" = ±oo. Deforming the contour of integra- 
tion to include a singularity, e.g. in the upper-half plane, 
we find that the integral has two contributions: one is 
provided by t" = 1 point, and another one is given by 
a combination of two paths that run from the complex 
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infinity to t" — 1 and backwards with the corresponding 
phase shift (Fig. [I}. As a result, we obtain: 



a+(oo) 



■ exp 



1 is 

7T 



2(5 
dx 
sinh : 



(27) 



exp[ 



• sinh(2x] 



where we changed the variables £5/ff = i cosh a; to sim- 
plify the expressions. In the limit S <C g 2 , the integral in 
the brackets is a constant equal to tr/6, which leads to 
the following asymptotic behavior of the transition prob- 
ability: 



|a+(oo)| 2 



9 



exp 



S 



(28) 



This expression correctly reproduces the exponential de- 
pendence of the transition probability on the Landau- 
Zener parameter g 2 /S. However, the exponential prefac- 
tor is larger than the exact value of unity [53| . The reason 
for this discrepancy is that the exponential dependence is 
a result of a delicate interference of the transition ampli- 
tudes in time, which can not be obtained within pertur- 
bative approach. Conversely, as we argue below, in the 
case of finite U or tf, the asymptotical behavior of the 
transition probability with S is analytic, allowing for a 
systematic treatment within the adiabatic perturbation 
theory. 

Let us now turn to such situation with finite ti or tt 
or both (positive or negative) and sufficiently small 5 . 
Then using Eq. ([20]) and ignoring the fast oscillating term 
we find 



«+(*/) f 



9 



9 



16 5 4 V (g 2 + A?) 3 (g 2 + X}) 3 



(29) 



In the case of A^ = —00 and Xf = 0, i.e. ti — —00 and 
tf = Q, this gives 



I6g 4 



(30) 



as it can be found solving exactly the LZ-problem (see 
Appendix IA")) . We would like to emphasize that in agree- 
ment with the general prediction (|20|) in the adiabatic 
limit the transition probability (|29p is quadratic in the 
rate S. We note that a more accurate asymptotic includes 
an additional exponential non-perturbative term (|24jl if 
A, and A/ have opposite signs With this addi- 

tional term Eq. (|2"9")) will give the correct asymptotic for 
\a + (tf)\ 2 , even when both |Aj| and |A/| are large. 



III. ADIABATIC DYNAMICS IN GAPLESS 
SYSTEMS WITH QUASIPARTICLE 
EXCITATIONS 

In the previous section we applied the adiabatic per- 
turbation theory to a two-level system with a time- 



dependent gap separating the eigenstates of the Hamil- 
tonian. Now we are interested in extending the analysis 
to the case of a many-particles system. If we consider a 
situation in which the system (initially prepared in the 
ground state) is characterized by gapped quasiparticle 
excitations, then clearly the previous analysis applies to 
each mode independently. Then under slow quench dif- 
ferent quasiparticle states can be excited. Since the tran- 
sition probability to each quasiparticle state quadrati- 
cally depends on the quench rate 5, we can expect that 
the density of created quasiparticles and other thermo- 
dynamic quantities will also quadratically depend on S. 
The situation becomes more complicated and more inter- 
esting if we consider a system with gapless excitations. 
Then for the low energy states the adiabatic conditions 
are effectively always violated (unless we consider quench 
rates which go to zero with the system size) and in princi- 
ple the adiabatic perturbation theory and the quadratic 
scaling can break down. Although the adiabatic pertur- 
bation theory in this case is no longer quantitatively cor- 
rect, here we will show that it can nevertheless be very 
useful in finding the scaling of various quantities with the 
quench rate and when this scaling becomes quadratic. 
Let us assume that we are dealing with a homogeneous 
system of quasiparticles characterized by a dispersion re- 
lation: 



s k = c(A)fe 2 



(31) 



where z is the dynamical exponent and c(A) is a pre- 
factor depending on the external parameter A. We as- 
sume that A changes linearly in time, A = St, between the 
initial and the final value A^ and A / respectively. The pa- 
rameterization pip explicitly demonstrates that during 
the evolution the details of the quasiparticle spectrum do 
not change, i.e. the system remains gapless and the ex- 
ponent z stays the same. Another interesting possibility 
where the system crosses a singularity such as a critical 
point, which violates this assumption, will be considered 
in the next section. In this work we also limit ourselves 
to global uniform quenches, where the coupling A is spa- 
tially independent. 

In this section we consider the situation where the 
ground state corresponds to the state with no quasi- 
particles. Such situation naturally appears in a variety 
of physical systems, e.g. bosonic systems with short- 
range interactions, fcrmionic systems without Fermi sur- 
faces such as gapless semiconductors (55|, graphene [56[ 
at zero voltage bias, some one-dimensional spin chains, 
that can be mapped to systems of fermions using the 
Jordan- Wigner transformation [57], Hfl , system of hard- 
core bosons in one-dimension (or Tonks gas) [H, [H, [6(| 
and so on. The situation with Fermi systems with the 
ground state corresponding to the filled Fermi sea re- 
quires special attention and will not be considered here. 
We do not expect any qualitative differences in the re- 
sponse to slow quenches in these Fermi systems. 

Dimensional analysis allows us to estimate the scal- 
ing of various thermodynamic quantities with the quench 



6 



rate. On the one hand, in a gapless system the transi- 
tions to the low energy states are unavoidable since any 
change in the coupling A looks fast (diabatic) with re- 
spect to them. On the other hand transitions to the high 
energy states are suppressed for small S because of the 
fast-oscillating phase factor entering the expression for 
the transition amplitude (I16p , so that the dynamics with 
respect to these states is adiabatic. It is straightforward 
to estimate the boundary separating diabatic and adia- 
batic states. Namely, if throughout the evolution £fc(A) 
becomes larger or comparable to e|(A), then the corre- 
sponding energy level is diabatic and the quasiparticles 
are easily created. On the other hand if £fe(A) -C £fc(A) 
in the whole interval of A € [Aj, A/] then the transitions 
are suppressed. Of course the implicit assumption here is 
that there is no kinematic constraint preventing creation 
of quasiparticles, i.e. that the matrix element for the 
transition is non-zero. In spatially uniform systems sin- 
gle quasiparticles typically can not be created because of 
the momentum conservation, so that quasiparticles can 
be created only in pairs with opposite momenta. This im- 
plies that one should use 2efc(A) instead of £fc(A) in the 
argument, however, this extra factor of two is not impor- 
tant for a qualitative discussion. The crossover energy 
£fe(A) is found by equating the two expressions: 



k ax 6 



c(A) 2 P 



(32) 



or alternatively 



£fc(A) ~ j dlnc(A) _ i QjneW 



ox 



c(A) dX 



More precisely one needs to find the maximal momen- 
tum satisfying the equation above in the whole interval 
of A e [Xi,Xfj. If there are no singularities during the 
evolution the derivative 9,\lnc(A) remains non-zero and 
bounded. Thus we see that the crossover energy £fc(A) 
scales linearly with the rate S. Then we can estimate, 
for example, the total number of quasiparticles created 
during the process as 



p(e)de 



d d k 



k<k 



(27T) 



d/z 



(34) 



We expect this scaling to be valid only when d/z < 2. 
Otherwise (as we will argue below) the contribution from 
high energy quasiparticles with e > e will dominate, re- 
sulting in the quadratic scaling of n ex . 

The result (|34|) can be derived more accurately using 
the adiabatic perturbation theory. Namely, let us per- 
form the scaling analysis of Eq. (JTSJ) ■ Because quasipar- 
ticles can be typically created only in pairs with opposite 
momenta one should use twice the quasiparticle energy 
in the dynamical phase O n in Eq. (|16p . The number 
of quasiparticles should be also multiplied by two, but at 
the same time the sum over momenta should only go over 
half the available states to avoid double counting of pairs 



with momenta k and —k. So this second factor of two 
can be absorbed extending the integration over momenta 
to the whole spectrum. Thus we get 



l r , , M|2 f d d k 



dX (£;|d A |0)exp 



2ik z 



c(X')dX' 



Ai 



(35) 



Rescaling the momentum k as k — 8 x l z r\ we find 



d d rj 
(2tt) c 



dX (i]6 1/z \dx\0}e 



2i v * / A A c(A')dA' 



If the integral over rj converges at large r\ and the matrix 
element (fc|9x|0) goes to a non-zero constant as k — > 0, 
then we get the desired scaling. The second condition 
means that there is no kinematic suppression of the tran- 
sitions to the low energy and momentum states. The ab- 
sence of this suppression was implicitly assumed in the 
elementary derivation of Eq. (j3~4"|) . The first condition 
of convergence of the integral in Eq. (|36j) over r\ implies 
that only low energy modes contribute to the total num- 
ber of generated quasiparticles since k ~ T)b x l z . This 
happens only in the case where d/z < 2, otherwise the 
quadratic scaling n cx ~ 5 2 , which we derived in the previ- 
ous section, coming from excitations to all energy scales, 
is restored. Technically the crossover can be seen from 
the fact that the transition probability to the state with 
momentum k at large k 3> S 1 ^ corresponding to rj 3> 1 
scales as l/rj 2z . This follows from combining Eqs. (|20|) 
and ([3"T]) . Clearly then only if d < 2z the integral over 
rj in Eq. (|36|) converges at large rj. Therefore after the 
rescaling k = 5 1, /z "q the upper limit of integration over rj 
can be sent to infinity. Otherwise the quasiparticle exci- 
tations with large momenta independent of 8 dominate 
n ex and we obtain the quadratic scaling according to the 
general result (|20[) . 

One can check that under the same conditions the ex- 
cess energy or heat [61] per unit volume in the system 
scales as: 

*2 ~ Jd^2 £ k( x f)\ a k(^f)\ 2 = \S\^c(X f ) 



d d r] , 



(27r) d 



dX (^V'l^io). 



2irf c(A')<U' 



(37) 



This scaling is now valid provided that (d + z)jz < 2. 
Otherwise the energy absorbtion comes from the high 
energy states with k S 1 ^ and in this case Q oc S 2 . 

Since we assume that the quasiparticles are created 
independently in different channels, i.e. the probabil- 
ity to create a pair of quasiparticles with momentum k 
is uncorrelated with the probability to create a pair of 
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quasiparticles with momentum k' ^ k, we can easily find 
the scaling of the (diagonal) entropy density of the sys- 
tem ffti: 



(38) 



where the factor of 1/2 comes from the fact that we need 
to account each state characterized by momenta q and 
—q once and sum only over half of the momentum modes 
to avoid over-counting. As before one can instead sum 
over all momentum modes but multiply the result by 1/2. 
In the domain of validity of the adiabatic perturbation 
theory |afc(A/)| 2 -C 1, the expression for the entropy is 
very similar to the expression for the number of quasi- 
particles. The extra logarithm clearly does not affect the 
scaling with 6 and we thus expect that Sd oc \S\ d ^ z for 
d/z < 2. 

Example: harmonic chain. Let us consider a specific 
example from Ref. [3l|, where one slowly changes the 
mass of particles in a harmonic chain. Specifically we 
consider the following Hamiltonian 



n h 2 



+ —\Pk\ , 



(39) 



where Xk and pk are conjugate coordinate and momen- 
tum. We will assume that k, playing the role of an inverse 
mass, linearly increases in time: n(t) = + St. For sim- 
plicity we will also assume that k/ S> k^. In Ref. [3lT ] a 
more complicated situation with k dependent on k was 
analyzed, however this fc-dependence is only important if 
Ki — > corresponding to a singularity in the spectrum. In 
this section we are interested in non-singular situations, 
therefore we assume that Ki is finite and that an extra 
possible dependence of k on k is irrelevant. The problem 
was analyzed both exactly and perturbatively in Ref. [3l| . 
Here we briefly repeat the analysis and present additional 
results. Within the adiabatic perturbation theory chang- 
ing k in time generates pairs of quasiparticle excitations 
with opposite momenta. It is straightforward to check 
that the matrix element (k,— k\d K \0) = 1/(4v / 2k) is in- 
dependent of k satisfying the second requirement (see the 
sentence after Eq. (j3l)|) ) necessary to get the correct scal- 
ing of the excitations. Substituting this matrix element 
together with the dispersion ek = k^Jnp s into Eq. (fT6]) . 
we find that within the perturbation theory the proba- 
bility to excite a pair of quasiparticles with momenta k 
and —k is 



1 

32 



"'t dn 



■ exp 



35 



3/2 



1 

72 



4 kKi jQgKi 

r ( o,-i- Y 



(40) 



where T(0,x) is the incomplete T-function. In the equa- 
tion above we used the fact that Kf 3> Ki so the upper 
limit in the integration over k can be effectively extended 



to infinity. The expression depends on the single dimen- 
sionless parameter: 

Cfe ~ jfc? 



where 



When ^k 3> 1, which corresponds to the high momen- 
tum modes, we have 



i 

128 k 2 



k 2 



and in the opposite limit Eq. (|40|) gives 



(41) 



(42) 



Hence at high energies the transition probability is pro- 
portional to 5 2 , as expected from the discussion in the 
previous section, while at small momenta k — > the tran- 
sition probability diverges. This is of course unphysical 
and indicates the breakdown of the perturbation theory. 
Note, however, that because of the small prefactor 1/72 
this divergence occurs at very small values of 

This problem can also be solved exactly (see details in 
Ref. [3l|). The initial ground state wave function can be 
written as a product of Gaussians: 



=n 



i 



(2^(^))V4 



exp 



Ffc 



4,7° (Ki 



> (43) 



where 



One particularly useful property of Gaussian functions is 
that for arbitrary time-dependence of n(t) (or p s (t)) the 
wave function remains Gaussian with Ok satisfying the 
first order differential equation: 



. dak (t) 
dt 



2p 8 k 2 a%{t) - -K[t). 



(44) 



The solution of this equation which satisfies the proper 
initial condition can be expressed through the Airy func- 
tions. It is convenient to introduce the effective width of 
the wave function at time tf. 

1 



It can be shown that [3111 



i + IM 

23/ fe 



(45) 
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where 



fk 



*Bi'(-£ 



/^Ai^f 5 )-^-^ 3 )' 



(46) 



As it occurs in the perturbative treatment, the solution is 
expressed entirely through the single parameter It is 
easy to show that the average number of excited particles 
with momenta k and — k in the Gaussian state charac- 
terized by the width cr£ is 



1 

n k = - 



4 



(47) 



This expression has the following asymptotics: for ^> 1 



1 



"3 



and for < 1 



3 2/3 r 2 (1/3) ^ 



(48) 



(49) 



For large momenta the exact asymptotic clearly coincides 
with the approximate one pi"j) (noting that n k = 2\a k \ 2 ) 1 
while at low energies as we anticipated the adiabatic per- 
turbation theory fails. The perturbative result clearly 
underestimates the number of excitations. This fact is 
hardly surprising because the adiabatic perturbation the- 
ory neglects the tendency of bosonic excitations to bunch 
together, leading to the enhancement of the transitions. 
We can anticipate that the result will be opposite in the 
fermionic case and in the next section we will see that 
this is indeed the case. 

Density of quasiparticles. The density of quasiparticles 
created during the process can be obtained by summing 
rik over all momenta: 



d d k 



(2*0 



(50) 



We remind that the factor of 1/2 is inserted to avoid 
double counting of pairs with momenta k and —k. In low 
dimensions d < 2 this sum is dominated by low momenta 
k ~ ks and the upper limit can be send to oo: 



K 8 



(2tt) c 



MO- 



(51) 



This gives the correct scaling n cx ~ 5 d / z anticipated from 
the general argument since z = 1. On the other hand 
above two dimensions the integral over k is dominated 
by high momenta close to the high energy cutoff and we 
have 



1 
128 



d d k 1 
(27r) d ¥' 



(52) 



This result again confirms our expectations that the num- 
ber of excited quasiparticles scales as S 2 in high dimen- 
sions. Note that when d = 2 the integral above is still 



valid but it should be cutoff at small k ~ kg, leading 
to the additional log-dependence of n cx on 8. We point 
out that the perturbative analysis predicts a very simi- 
lar qualitative picture: it correctly predicts the density 
of excitations in high dimensions d > 2 and gives the 
correct scaling in low dimensions d < 2. However, in 
the latter case the perturbative analysis gives a wrong 
prefactor. We note that there are situations when the 
adiabatic perturbation theory can fail completely pre- 
dicting even incorrect scaling (3l| . This can happen, for 
example, if the initial coupling k is very small. Then 
the integral ([50]) can become infrared divergent at small 
momenta and should be cutoff by the inverse system size 
or another large spatial scale (e.g. the mean free path). 

The smallest (and the only) physical dimension where 
the scaling n ox cx \S\ d ^ z is valid in our situation is d= 1. 
Accurate evaluation of the quasiparticle density in this 
case gives 



w 0.0115^. 



(53) 



For completeness we also quote the perturbative result 
obtained by integrating Eq. (|4"0"]) : 



port 



omoiks. 



(54) 



Clearly the difference between the exact and perturbative 
results is very minor in this case. 

In two dimensions both perturbative and exact treat- 
ments give the same result 



1 



256tt 



k 2 \n 



ks 



(55) 



where A is the short distance cutoff. At higher dimen- 
sions, as we mentioned above, results of perturbative and 
exact treatments are identical and non-universal: 



CkU d ~ 2 . 



(56) 



It is interesting to note that adiabatic conditions are 
determined now by the high momentum cutoff: S ~ 
1/A d / 2_1 . This sensitivity to A might be an artifact of 
a sudden change in the rate of change of the parameter 
k resulting in the infinite acceleration A at initial (and 
final) times. However, analyzing this issue in detail is 
beyond the scope of this work. 

Heat (excess energy). The quasiparticle density is not 
always an easily detectable quantity, especially if the sys- 
tem is interacting and the number of quasiparticles is not 
conserved. A more physical quantity is the energy change 
in the system during the process, which is equal to the 
external work required to change the coupling n. This 
energy consists of two parts: adiabatic, which is related 
to the dependence of the ground state energy on k, and 
heat, i.e. the additional energy pumped into the system 
due to excitations of higher energy levels. The first adi- 
abatic contribution depends only on the initial and final 
couplings Ki and k/ but not on the details of the process. 
On the contrary the heat (per unit volume) Q is directly 
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related to the rate 8. The microscopic expression for Q 
can be obtained by a simple generalization of Eq. (I5C 
as shown in Eq. ([37| : 



Q 



d d k 
(2n) d 



■knk- 



(57) 



This integral converges at large k only for d < 1. There- 
fore in all physical dimensions it is dominated by the 
high energy asymptotic (|4"Tj) so Q oc S 2 . In one dimen- 
sion there is an extra logarithmic correction 



(58) 



Entropy. Finally let us evaluate the generated diagonal 
entropy in the system. The latter is formally defined as 
(see Ref. [H) 



Sd 



7? 5> 



log pr, 



(59) 



where p nn are the diagonal elements of the density matrix 
in the eigenbasis of the (final) Hamiltonian. At finite 
temperatures it is this entropy which is connected to heat 
via the standard thermodynamic relation AQ — TASd- 
However, at zero initial temperature this relation breaks 
down and one should analyze Sd = ASd separately. 

Since in our problem different momentum states are 
decoupled the d-entropy is additive: 



Sd 



—y 

2L d ^ 



1 

Sk = 2 



d d k 
(27r) d 



Sk, 



(60) 



where the factor of 1/2 is again present in order to 
avoid double counting of quasiparticle excitations cre- 
ated in pairs. Within the adiabatic perturbation theory 
only one pair of quasiparticles can be excited. Therefore 
s P Glt ~ _ \ak\ 2 In | CKfc | 2 - Note that at large k the entropy 
per mode falls off with k essentially in the same manner 
as \otk\ 2 so we conclude that the entropy is dominated by 
small energies below two dimensions. Then we find 



gpert 



1 



d d £, 



HOI 2 m KOI 2 



(61) 



In dimensions higher than two the entropy is dominated 
by high momenta so 



s 



pert 



i 



d d k 1 



In 



k'i 



256'"° J {2ix) d k 2 
In one dimension Eq. (|6ip gives 

SSff « 0.02^ 



Vl28fc s 



(62) 



(63) 



To calculate the exact value of the entropy density one 
needs to project the Gaussian wave function describing 
each momentum mode to the eigenbasis and calculate the 



sum (|59|) . There is, however, a simple shortcut allow- 
ing to instantly write the answer. One can easily check 
that the probabilities to occupy different eigenstates are 
identical to those of the equilibrium thermal ensemble. 
Therefore the entropy per mode can be expressed through 
the average number of excited quasiparticle pairs n k /2 as 

s k = -(n fc /2) ln(n fc /2) + (1 + n k /2) ln(l + n k /2). (64) 

Using explicit solution for n k we find 

S d =i « 0.026fc 5 . (65) 

The result is again quite close to the perturbative one. 

In two dimensions the entropy density is also readily 
available from the expressions above: 



Sd=2 



5127T 



(66) 



IV. 



ADIABATIC DYNAMICS NEAR A 
QUANTUM CRITICAL POINT 



Let us apply the analysis of the previous section to the 
case of crossing a quantum critical point. As before we 
consider the situation in which the system is prepared in 
the ground state, characterized by some initial coupling 
Xi, then this coupling is linearly tuned in time until a 
finite value A/. We assume that the system undergoes 
a second order quantum phase transition at A = 0. We 
differentiate the two situations: (i) when Ai is finite and 
negative and Xf is finite and positive and (ii) when cither 
Ai = and A/ is large and positive, or A/ = and Xi 
is large and negative. As in the previous section we will 
first give a general discussion closely following Ref. [27j 
and then analyze specific examples. 



A. Scaling analysis 

The non-adiabatic effects are especially pronounced in 
the vicinity of a quantum critical point (QCP), where 
one can expect universality in the transition rates. Near 
a QCP for A ^ the system develops a characteristic 
energy scale A, vanishing at the critical point as A ~ 
|A| zl/ , where z is the dynamical exponent, and v is the 
critical exponent of the correlation length, £ ~ \X\~ V [58|. 
This energy scale can be either a gap or some crossover 
scale at which the energy spectrum qualitatively changes. 

As before, let us perform the scaling analysis of 
Eq. (|16p assuming that the quasiparticles are created in 
pairs of opposite momenta. First we rewrite the dynami- 
cal phase factor entering the expression for the transition 
amplitude (HHJ) as 



e fe (A)-e (A) 



dX'(e k (X')-e a (X')). (67) 
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Near the QCP the quasiparticle energy can be rewritten 
using the scaling function F as: 



e k (X)-e (X)=X zu F(k/X v ). 



(68) 



For x ^> 1 the function F(x) should have asymptotic 
F(x) oc x z , reflecting the fact that at large momenta 
the energy spectrum should be insensitive to A. At small 
x, corresponding to small k, the asymptotical behavior of 
the scaling function F(x) depends on whether the system 
away from the singularity is gapped, then F(x) — >const 
at x — > 0, or gapless, then F(x) oc x a with some positive 
exponent a. The scaling ([68]) . inserted in (|67|) . suggests 
the change of variables 



(69) 



In addition to analyze Eq. JI6 
ansatz for the matrix element 



we adopt the scaling 



{k\dx\0) 



(k\V\0) 



e*(A)-E (A) A 



G(k/xn, 



(70) 



where G(x) is another scaling function. This scaling di- 
mension of the matrix element, which is the same as the 
engineering dimension 1/A, follows from the fact that 
the ratio of the two energies (fc|AV|0) and £fc(A) — £o(A) 
should be a dimensionless quantity. At large momenta 
k X" this matrix element should be independent on A 
so G(x) oc x' 1 /" at x ^> 1. This statement must be true 
as long as the matrix element (fc|V|0) is non-zero at the 
critical point, which is typically the case. In the opposite 
limit x -C 1 we expect that G(x) oc x@ , where (3 is some 
non-negative number. 

These two scaling assumptions allow one to make some 
general conclusions on the behavior of the density of 
quasiparticles and other thermodynamic quantities with 
the quench rate for the adiabatic passage through a quan- 
tum critical point. Thus 



d d k 
{2n) d 



(71) 



where, after the rescaling ([69 



a(rj) 



dt±G ( f 



exp 



(72) 

Note that if Ai < and A/ > then £j is large and nega- 
tive and £/ is large and positive. If we start (end) exactly 
at the critical point then £j =0 — 0). Note that the 
integral over £ is always convergent because: at large £ we 
are dealing with a fast oscillating function, and at £ ~ 
there are no singularities because of the scaling properties 
of G(x). The only issue which has to be checked carefully 
is the convergence of the integral over rj at large r\. If this 
integral converges then Eq. (|TTj) gives the correct scaling 
relation for the density of quasiparticles with the rate S 
found in earlier works [1?], HI] . If the integral does not 



converge then the density of created defects is dominated 
by high energies and from general arguments of Sec. [TT] 
we expect that n cx oc S 2 . As it is evident from Eq. (|71[) . 
the crossover between these two regimes happens when 
d = 2z + 2/v [7l|. To see how this crossover emerges from 
Eq. (|71j) we analyze the asymptotical behavior of a(rj) at 
77 3> £ v . In this limit the integral over £ can be evaluated 
in accord to the discussion given in Sec. [TT] because the 
exponent in Eq. (|72|) is a rapidly oscillating function of 
£. Using the explicit asymptotics of the scaling functions 
F(x) and G(x) at large x we find 



a(rf) oc 



1 



V 



z+l/v ' 



(73) 



Then the integral over 77 in Eq. (|TTj) converges at large 77 
precisely when d < 2z + 2 /v. Instead when the integral 
does not converges, the scaling n ex oc \5\ dv / ( zvJrl ) breaks 
down and we get n ox oc S 2 . As in the previous section 
we can expect logarithmic corrections at the crossover 
between these two scaling behaviors. 

One can similarly analyze the dependence of the heat 
Q and the diagonal entropy S4 on the rate 6. The entropy 
has the same scaling as the number of quasiparticles. We 
note that in general the heat is universal only if the pro- 
cess ends at the critical point A/ = 0, where e k oc k z . 
Then it is easy to see that 



<x\8\- 



(74) 



This scaling is valid for d < z + 2/v, while it becomes 
quadratic when d > z + 2jv. If the final coupling A/ is 
away from the critical point then the ^-dependence of Q 
becomes sensitive to the behavior of the spectrum. Thus 
if there is a gap in the spectrum, Q has the same scaling 
as n cx and Sd, since each excitation roughly carries the 
same energy equal to the gap. If the spectrum at Ay is 
gapless then the scaling (|T4"1) remains valid. 



B. Examples 

Here we will consider several specific models illustrat- 
ing the general predictions above. In particular we will 
analyze the transverse field Ising model, which serves as a 
canonical example of quantum phase transitions (EH ] and 
which was used as an original example where the gen- 
eral scaling (fTTj) was tested [13, HH, ■ The transverse 
field Ising model also maps to the problem of loading one- 
dimensional hard-core bosons or non-interacting fermions 
into a commensurate optical lattice potential [351 ] and de- 
scribes the so-called Toulouse point in the sine-Gordon 
model, where this model can be mapped to free spinless 
fermions 1571. 
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1. Transverse field Ising and related models 

The transverse-field Ising model is described by the 
following Hamiltonian: 



Hi 



3 



(75) 



For simplicity we will focus only on the domain of non- 
negative values of the transverse field g. This model 
undergoes a quantum phase transition at g = 1 [58| 
with g > 1 corresponding to the transversely magnetized 
phase and g < 1 corresponding to the phase with lon- 
gitudinal magnetization. It is thus convenient to use 
A(i) = g(t) — 1 as the tuning parameter. Under the 
Jordan- Wigner transformation: 

a] = -( Cj + c]) H(l - 24a), o? = l- 2c]c 3 (76) 

i<3 

the Hamiltonian assumes the quadratic form: 

Hi = - X] c ] c i+! + 4+i c J + c l c i+i + c J+ lC J ~ 2 9(t)4 c 3 

(77) 

and can be diagonalized using the Bogoliubov transfor- 
mation in the momentum space: 



Cfe = Ik cos(6» fc /2) + i sm(9 k /2)jl k , 



where 



tan 8i. = 



sin(fc) 



cos(fc) — g(t) 

After this transformation the Hamiltonian becomes 



(78) 



(79) 



where 



£fc = 2yl + .g 2 - 2gcos(fc). 



(80) 



51) 



The ground state of this Hamiltonian factorizes into the 
product: 

|fio) = II (coB(d fc /2)+iBin(flfc/2)4cL fc ) |0). (82) 



The excited states can be obtained by applying various 
combinations of operators 7^ to the ground state above. 
However, because of the momentum conservation, only 
the excited states obtained by acting on the ground state 
by the products T k T_ k are relevant. Because the exci- 
tations to different momentum states are independent, 
the problem effectively s plit s into a sum of independent 
Landau-Zener problems [29] and can be exactly solved. 
In the case when gi and gf lie on the opposite sides of 



the quantum critical point the transition probability in 
the slow limit is approximately given by [29j 



Pk ~ exp 



2nk 2 



(83) 



The density of the excited quasiparticles is then 



1 
2H 



Pkdk 



2ttV2 



O.llVd. (84) 



Similarly one can find the entropy density generated dur- 
ing the process 

1 r°° 

Sd~-~r I dk \pklnpk + (1 -pfc)ln(l -pk)\ 
w 0.052V5. (85) 

Both the expressions for n GX and Sd agree with the gen- 
eral scaling law (|7"Tj) with d = z = v = 1. Since in 
the final state all the excitations are gapped the heat in 
this case is approximately equal to the number of excited 
quasiparticles multiplied by the gap in the final state and 
thus has the same scaling as n cx . 

The same problem can be also solved using the adia- 
batic perturbation theory. Let us note that as it follows 
from Eq. (|52"|) the transition matrix element reads [27j 

sink 



(k,-k\dx\0) = (k,-k\d g \0) = - 



2 1 + g 2 - 2gcos(k)' 

, (86) 

In the limit of small S only the transitions happening in 
the vicinity of the critical point contribute to n ox . In this 
case we have 



{k,-k\d x \0) 



(87) 



2 A 2 + k 2 ' 

which clearly satisfies the scaling (|70l) . Under the same 
approximation we can use that Sk ~ 2\Jk 2 + A 2 , which 
in turn satisfies the energy scaling (|68|) . Substituting 
these expansions into Eq. (|16p and extending the limits 
of integration over A to (—00, 00) we find that 



Oik 



d\- 



k 



A 2 + k' 



■ exp 



V / dX'Vk 2 
5 Jo 



A' 2 



where the additional factor of two in the exponent comes 
from the fact that two quasiparticles are created during 
each transition. It is now straightforward to evaluate the 
perturbative expression for n^* . 



pert 



dk, 
2^ 



2 M 0.21VS. 



(89) 



Note that unlike the bosonic case discussed in the previ- 
ous section, the adiabatic perturbation theory now over- 
estimates the number of created quasiparticles. This hap- 
pens because this perturbation theory does not take into 
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account the Pauli blocking which prevents more than one 
pair of quasiparticles with momenta fc, — fc to be excited. 
Technically the perturbative transition probability |afc| 2 
can exceed unity. This prevents from computing the gen- 
erated entropy using Eq. (f38|) because there will be spu- 
rious negative contributions. To get a sensible expression 
one needs to integrate only over the momenta satisfying 
I cifc 1 2 < 1, and this results in: 



-fpcrt 



1 

2^ 



dk\a k \ 2 In |a fc | 2 



0.022v^. (90) 



The perturbative argument can be somewhat improved 
explicitly using the fact that the quasiparticles are 
fermions and adding an additional contribution com- 



ing from (1 



iln(l — \au\ )■ In this case S% 



pert 



0.03&VS, which is closer to the exact result. 

In a similar spirit we can consider the situation when 
either A; = or A/ = 0, i.e. when the initial or final state 
exactly corresponds to the quantum critical point. This 
situation requires fine tuning from the point of view of 
crossing a quantum phase transition. However, it natu- 
rally appears in other contexts. For example, the prob- 
lem of loading hard core bosons into a commensurate 
periodic potential exactly describes this situation [35| . 
Since we consider A E (0, oo) then Eq. ([55]) becomes: 



Oik 



dX 



A 2 + fc 2 



exp 



(91) 

then it is easy to see that nPx rt (0,oo) = ^ng£ rt (— oo, oo), 
therefore the scaling of the density of quasiparticles (and 
entropy) remains the same n ex cx y/5. This problem can 
be also solved exactly and the results remain very close to 
the perturbative case. We will postpone a careful analysis 
of this problem to the next section. 



2. sine-Gordon model: Toulouse point and limit of free 
massive bosons 



A convenient playground to test the general scaling 
laws presented in Sec. IIV Al is the sine-Gordon model 
(SG) - a one dimensional model described by the Hamil- 
tonian: 

H = i J d x [U{x) 2 + (d x (t)) 2 - 4A cos(/3 0)] . (92) 

Here H(x) and <j){x) are conjugate fields, A is the tun- 
ing parameter and (3 — 2\/ttKsg is a constant. From the 
renormalization group (RG) analysis it is known that the 
cosine-term is a relevant perturbation to the quadratic 
model only if < Ksg < 2, and therefore the system is 
gapped at any finite A [57|], while for Kg G > 2 the sys- 
tem remains gapless. Turning on the interaction in the 
regime -Ksg < 2 is akin starting at the critical point and 



driving the system into the new gapped phase. Hence the 
scaling of the density of excitations and other quantities 
should be described by the critical exponents in accord 
to Eq. (fTTj) as it was indeed shown in Ref . [35| . The spec- 
trum of the SG model consists of solitons and antisolitons 
for 1 < Ksg < 2 and in addition to those for Ksg < 1 
there are also breathers excitations. The point Ksg = 1 
is the Toulouse point described by free fermions and in 
the limit Ksg — > the SG model effectively describes a 
system of free massive bosons with B\ breathers being 
the only surviving excitations (see e.g. Ref. 63]). The 
energy spectrum for each momentum fc is of the form: 



£fc = V'fc 2 + ' 



(93) 



where m s is the soliton (or breather) mass, that scales 
with the external parameter as m s ~ A 1 ' ( 2 -- K "so)_ There- 
fore it is evident how this model fulfills the assump- 
tion (|68[) with the critical exponents z = 1 and v — 
1/(2- K SG ). 

The SG model gives the correct low energy description 
of: (i) interacting bosons in a commensurate periodic 
potential and (ii) two one-dimensional condensates (Lut- 
tinger liquids) coupled by a tunneling term [35|. In the 
former case increasing A in time corresponds to loading 
bosons into an optical lattice and in the latter case in- 
creasing A describes turning on the tunneling coupling. 
We note that <fi describes the density modulation in the 
first situation and the relative-phase between the two su- 
perfluids in the second. The SG model also natu rally 
appears in several other one-dimensional systems [5 71 ]. 
To simplify the analysis we will consider only the situa- 
tion where either A^ = and Xf large or vice versa. For 
a generic value of Ksg the problem can only be solved 
perturbatively [35j and the corresponding discussion is 
beyond the scope of this work. Here we will discuss only 
two specific solvable limits Ksg <C 1 and Ksg = 1. 

a. Ksg ^ 1 ; / r ee bosons. This limit can be nat- 
urally realized in the situation of merging two weakly 
interacting one-dimensional condensates. Then Ksg = 
1 / (2K) , where K » 1 is the Luttinger liquid parameter 
describing the individual condensates [35j ]. Since Ksg 
and hence (3 in the Hamiltonian (|92[) is small we can ex- 
pand the cosine term and get a quadratic Hamiltonian 
that in the Fourier space has the form: 



H = 



Kk{t)\<t>k\ 



(94) 



with Kk(t) = fc 2 + 2A(i)/3 2 , where as before we assume 
that X(t) — St. We see that in this case the problem is 
equivalent to the harmonic chain ([5^1) already considered 
in Section IIIII T he on ly difference is in the excitations 
spectrum e k = \j Kk(t) which is now gapped: 



efc = \Jk 2 + cA, 



(95) 



with c = 2/3 2 . The perturbative and exact solutions fol- 
low the similar steps showed in Sec. IIIII In the pertur- 
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bative case we find 



1 

32 



K i dn 



■ exp 



4i 

3c$' 



,3/2 



1 

72 



4fc 3 



where we used the fact that Kj = fc 2 . 
introduce the rescaled momentum 

fc«5 



where 



fc* = (c^) 1 / 3 . 



(96) 

It is convenient to 
(97) 

(98) 



Then all the results can be found from those of the har- 
monic chain (cf. Eq. (|40p ), performing the mapping 

6 < — > Cf- 

In the limit ^ 1, corresponding to the high momen- 
tum modes, the transition probability becomes 



1 



128C 6 



1 (cS) 2 
128 "P" 



and in the opposite limit 



1 

72 



^ I In Cf I 



(99) 



(100) 



For the exact solution of the problem we again follow 
the Gaussian functions ansatz as in Eq. (|43jl , with initial 
value 



V«*(0) 



and time dependence satisfying the equation: 

.^(t) =2 ^ 1 
eft few 2 w 



(101) 



The solution of this equation is analogous to the one of 
Eq. (14~4"1) with the only difference that the function //; (cf. 
Eq. (14"6"1) ~) now becomes 



/fc 



qBi(-Cg)-^(-Cg) 

C*Ai(-C2)-iAi'(-C2) : 



(102) 



which gives the following asymptotics for the average 
number of excited particle pairs with momenta fc and 
-fc: 



1 



1 {cSf 

64C| ~ 64 k 6 



for Cfc ^ 1 and 



(cS) 



1/3 



3 2 / 3 r 2 (i/3)a 3 2 / 3 r 2 (i/3) 



(103) 



(104) 



in the opposite limit. 

Density of quasiparticles. We note that the exact re- 
sult gives weak logarithmic divergence of the density of 
quasiparticles with the system size coming from the 1/fc 
dependence of n k at small k (see Eq. (|104[) ): 

nS act M (c^logfoL) w 0<03 3 (off)1/ai (fc L) 
2^3 2 r 2 (l/3) 

(105) 

while the perturbative solution gives: 



<« rt ra 0.068(c<5) 1/3 . 



(106) 



As previously found, the perturbative result under- 
estimates the number of excitated quasiparticles, since 
it does not take into account the bosonic enhancement 
of the transitions. Emergence of the length dependence 
in the expression for n cx indicates the breakdown of the 
adiabatic perturbation approach and corresponds to a 
different (non-adiabatic) response of the system accord- 
ing to the classification of Ref. [3l|. Physically this di- 
vergence comes from the effect of bunching of bosons at 
small momenta and overpopulation of low momentum 
modes, which can not be captured in the lowest order of 
the adiabatic perturbation theory. A similar analysis can 
be also performed for the heat and the entropy. In the 
situation in which the system starts at the critical point 
(Aj = 0) and A/ is finite, both the entropy and the heat 
have a very similar behavior as n GX , i.e. proportional to 
5 1 / 3 and showing a weak logarithmic divergence with the 
system size. In the opposite case where Xi is finite but 
A/ = the expressions for the density of excitations and 
entropy do not change, while the expression for the heat 
becomes different because each mode now carries energy 
£fc ~ k proportional to the momentum. This removes the 
logarithmic divergence and both perturbative and exact 
results give Q ~ fc 2 ~ (cS) 2/3 

b. Ksg = 1- Tonks- Girardeau gas. In the limit of 
Ksg — 1 the repulsive interaction between bosons is in- 
finitely strong, therefore the particles behave as impene- 
trable spheres (hard-core bosons). It is well known that 
in this limit the system known as the Tonks-Girardeau 
gas [HI can be mapped into an equivalent system of 
free spinless fermions (the corresponding limit of the SG 
model describes the so called Toulouse point). There- 
fore the dynamical problem of loading hard-core bosons 
into a commensurate optical lattice can instead be ap- 
proached with the much simpler analysis of free fermions 
in a periodic potential. 

To understand the dynamics in this case we need to 
solve the Schrodinger equation of free fermions in a peri- 
odic potential with time dependent amplitude V(x, t) = 
V(t) cos(2kfx), where fc/ = tt /a is the Fermi momentum 
and a is the lattice spacing. The potential V(t) is re- 
lated to the coupling A(t) in the SG Hamiltonian (|9"2"|) 
according to V(t)n = 4A(t), where n — kf/ir is the elec- 
tron density [H, [64|- We assume V(t) = 8yt, therefore 
the rate 8y is related to the rate 5 by simple rescaling: 
5y = 45 /n. We restrict the analysis to the two lowest 
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bands of the Brillouin zone (Fig. [2]) and we linearize the 
spectrum close to the Fermi momentum. These two ap- 
proximations, justified for analyzing the low energy exci- 
tations that we are interested in, make the fermion prob- 
lem equivalent to the SG problem. It follows that for each 
momentum k the problem is described by a Landau-Zcncr 
Hamiltonian 



H k = 



V(t)/2 A fe 
A, -V(t)/2 



(107) 



where = (e& — ek-2kt)/2 is half the energy difference 
between the 1st and the 2nd band. In the linearized 
approximation (under which the mapping to SG model 
is valid) we have ek ~ v f (k — kf), where Vf is the Fermi 
velocity. Choosing the units where Vf = 1 makes the free 
fermion problem identical to the SG problem with the 
Hamiltonian 



V 


w 







deep lattice 
shallow lattice 



o o o 
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To find the total density of excited particles we have to 
sum the probability p ex (k) over all the momenta in the 
first Brillouin zone: 



£ Pex(fc). 
fe£[— 7r/a,7r/a] 



(110) 



It is convenient to introduce a shifted momentum q = 
(k — kf). Then we have Tfe = (q/kg) 2 with kg — *^L, In 
the limit of small 5y the upper limit of the integral over 
q can be sent to infinity, therefore 



dq 
2^ 



PeAq 2 /k 2 5 ) » 0.12fe. 



(Ill) 



We see that the exact solution confirms the general scal- 
ing n ex ex y/S (see Eq. (fTTj) ) provided that d = z = v = 1. 

This problem can be also solved perturbatively. In fact 
the loading problem becomes equivalent to the transverse 
field Ising model. For example from the mapping to the 
LZ problem it easy to see that the matrix element ap- 
pearing in Eq. (fT6|) is: 



(112) 



-k\d\\-k) = r- 



2 (2A fe ) 



V'2 



which is similar to Eq. (|87[) . The situation is also anal- 
ogous for the energy spectrum. The net result of the 
calculation is 



pert 



O.Ukg. 



(113) 



The results for the entropy and heat for the loading prob- 
lem have identical scaling with the loading rate. 



FIS. 2: TG-gas: the loading into a commensurate lattice 
problem can be mapped into a two-level system composed of 
lower filled band and upper empty band, the excitations are 
the particles hopping on the upper band (figure taken from 

HI). 

The problem of turning on the potential amplitude 
from zero maps then to a direct sum of half LZ-problem 
(f € [0,+oc)) (see Appendix IX)) and can be solved ex- 
actly (see also Ref. [65j]). Then the probability of exciting 
a particle from the lower to the upper band during the 
loading process is: 



Pox(Tfc) =1-2 



-*7VB _ fJTTk 

sinn I 

7TTfe V 4 



with 



Tfe 



4A| 
Sv ^ 



(108) 



(109) 



SUDDEN QUENCHES NEAR QUANTUM 
CRITICAL POINTS 



So far we focused exclusively on the adiabatic regime 
of small d. For simplicity we restrict the discussion in this 
section only to the situation where we start at the critical 
point and increase A in time: X(t) = St. We argued that 
in low dimensions, when dvj(zv + 1) < 2, the dynamics 
is dominated by the low energy excitations with energies 
e < e* = \5\ zp '< '( zu+1 ) . The energy scale e* corresponds 
to the value of the tuning parameter A* = | ( 5| 1 /( ziy + 1 ) 
(because of the general scaling relation e ~ A ZI/ ). So 
equivalently we can say that in low dimensions the non- 
adiabatic effects are dominated by the transitions occur- 
ring in the vicinity of the critical point |A| < A*. A sim- 
ple qualitative way to understand the adiabatic dynamics 
in this case is to split the time evolution into two do- 
mains: t < t* = X*/S - i/|<S|W(^+i) and t > t*. Then 
in the first domain the dynamics can be thought as ap- 
proximately fast (sudden) and in the second domain the 
dynamics is adiabatic (see also discussion in Ref. l35]). 
Then we can think about the slow quench as a sudden 
quench with the amplitude of the quench being equal 
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to A*. I.e. one can expect that the scaling of the den- 
sity of excitations and other thermodynamic quantities 
can be approximately obtained by projecting the initial 
ground state corresponding to the critical point (A = 0) 
to the eigenstates of the new quenched Hamiltonian with 
A = A*. Such approach was indeed successfully applied to 
the problem of quenching the system through the BCS- 
BEC crossover [66| . This argument implies that for sud- 
den quenches we should have that n cx ~ \X*\ dv (45[. For 
dv < 2 this is a non-analytic function of the quench rate, 
which can not be obtained within the ordinary perturba- 
tion theory. As we will show below this scaling, however, 
immediately follows from the adiabatic perturbation the- 
ory applied to sudden quenches. Also we will show that 
when dv > 2 this scaling fails and the perturbative ana- 
lytic result is restored n ex ~ (A*) 2 . As in the case of slow 
quenches the analytic quadratic dependence comes from 
the dominance of the high energy excitations [45|] . 

Let us return to the discussion of Sec. HH where we in- 
troduced the adiabatic perturbation theory. We empha- 
size that the word adiabatic only means that we are work- 
ing in the instantaneous (adiabatic) basis. The small pa- 
rameter in this theory is the probability to excite higher 
energy levels. For slow quenches the excitation proba- 
bility is small because the rate 5 is small, while for fast 
quenches this probability is small because the quench am- 
plitude A* is small. Thus for sudden quenches of small 
amplitude we can still use Eq. (fT6)) . with further simpli- 
fication that the phase factor O n — 8o — > since it is 
inversely proportional to the rate <5 — > oo. Thus instead 
of Eq. (|16p we can write [45| 



a n (A*)«- / dA'<n|fl A '|0)= / dX' 



i\V\0) 



E n (X>) - E (X>Y 
o o 

(114) 

Note that in the case where there is a finite gap in the 
spectrum, the difference E n {X') — Eq(X') remains large 
for all values of A' G [0, A*]. If the same is true for the 
matrix element then this expression reduces to the one 
from the ordinary perturbation theory: 



a„(A*) « A* 



(n\V\0) 



E n (0)-E o (0y 



(115) 



The advantage of using the expression qil4|) for a n over 
the standard perturbative result (| 1 1 5[) is that Eq. (|114p 
does not give an explicit preference to the initial state 
A = over the final state A = A* and it does not assume 
the analytic behavior of a n with A* . As we will see both 
these points are necessary for studying quench dynamics 
near QCP. 

Let us apply Eq. (|1 14[) to the problem of quenching 
starting from the critical point. Using the scaling ansatz 
for the matrix element (|70|) . we find that the density of 
excitations is approximately given by 



d d k 
{2ir) d 



dX 



G 



(116) 



To analyze this expression it is convenient to change vari- 
ables as: A = A*£, k — (X*) v r]. Then we obtain 



d d f] 

(2-K) d 



ld A G (j± 



(117) 



This expression gives the desired scaling n cx ~ (|A|*) dj/ 
provided that the integral over r\ converges at large r\. 
From the asymptotics of the scaling function G(x) oc 
at large x, discussed earlier, we see that the neces- 
sary condition for convergence of the integral is d < 2/v 
or dv < 2. In the opposite limit, dv > 2, the high energy 
quasiparticles dominate the total energy and one can use 
the ordinary perturbation theory (linear response) which 
gives 



*|2 



|A*| 



d d k |(fc|U|0) o | 2 
(2n) d \e a k -e^ 



(118) 



where all the quantities are evaluated at the critical point 
A = 0. We note that the quantity multiplying |A*| 2 is 
called the fidelity susceptibility [67| . Thus the regime of 
validity of the universal scaling (|117|) corresponds to a 
divergent fidelity susceptibility. One can make similar 
analysis of the scaling of the heat and the entropy: 



K | A *|(d+*V s d oc\X*\ dv . 



(119) 



As in the case of n ex these scaling laws are valid as long as 
the corresponding exponents are less than two. One can 
verify that these scalings are indeed reproduced for the 
models we analyzed in the previous section. Both adia- 
batic perturbation theory and the exact calculation give 
the same scaling, however, the perturbative approach 
gives a mistake in the prefactor. 

From the above analysis we see that in low dimensions 
dv < 2 there is indeed a direct analogy between slow 
and sudden quenches. One gets the same scaling if cor- 
rectly identifies the quench rate i5 in the former case and 
the quench amplitude A* in the latter: A* ~ |<5| 1 /( ZI/+1 ). 
This analogy has direct similarity with the original ar- 
guments by Kibble and Zurek who predicted the scaling 
(fTTj) in the context of topological defect formation while 
crossing classical phase transition with the parameter 5 
playing the role of temperature quenching rate [3 [4t| . 
The main argument in the Kibble-Zurek (KZ) mechanism 
is that there is a divergent relaxation time scale near 
the critical point and a corresponding divergent length 
scale. The topological excitations with distances big- 
ger than this length scale do not have a chance to ther- 
malize and remain in the system for a long time even 
if they are thermodynamically forbidden. A simple es- 
timate gives that this length scale behaves precisely as 
f - l/\5\ v /(* v+v > [H, resulting in a density of topologi- 
cal defects n cx ~ which is the same as in Eq. (fTTj) . 
One can thus think about the KZ mechanism as adiabatic 
quenching of the temperature up to the scale T*, corre- 
sponding to the correlation length £, followed by sudden 
quench of the temperature to zero, so that remaining 
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topological excitations essentially freeze. Essentially the 
same argument we used here in the quantum case and it 
indeed works qualitatively right for dv < 2. 

The situation is different in the regime dv / '(zv + 1) < 
2 < dv. In this case for slow quenches the universal 
scaling is still applicable, i.e. only the low energy exci- 
tations created at A < A* are important for the scaling 
of n cx while for sudden quenches this is no longer the 
case. Namely, the scaling of n ox with the quench am- 
plitude becomes quadratic and non-universal, i.e. sensi- 
tive to the high energy cutoff. Therefore in this regime 
the analogy between slow and sudden quenches becomes 
misleading. This makes an important difference with the 
Kibble-Zurek mechanism where such issues do not arise. 
If dv/(zv + 1) > 2 then the scaling becomes quadratic in 
both cases (with the rate 8 for slow quenches and with 
the amplitude A* for sudden quenches), with the main 
contribution to excitations coming from the high energy 
quasiparticles. 



VI. EFFECT OF THE QUASIPARTICLE 
STATISTICS IN THE FINITE TEMPERATURE 
QUENCHES 

One can try to extend the analysis we carried through 
(both for sudden and slow quenches) to the finite tem- 
perature situation. We note that because we consider 
isolated systems with no external bath, the temperature 
only describes the initial density matrix. To realize this 
situation one can imagine that the system weakly couples 
to a thermal bath and reaches some thermal equilibrium. 
Then, on the time scales of the dynamical processes we 
are interested in, this coupling has a negligible effect and 
the dynamics of the system is essentially Hamiltonian. 
One does not even have to assume a coupling to the ex- 
ternal bath if we arc dealing with ergodic systems, since 
they are believed to reach thermal equilibrium states by 
themselves. By now this problem remains unsolved in 
the most general case. As we will argue below at finite 
temperatures the statistics of the low energy quasipar- 
ticles plays a key role. In general one can have critical 
points where the low energy excitations have fractional 
statistics or do not have well defined statistics at all. 

In this section we will consider only a relatively sim- 
ple situation in which the statistics of quasiparticles is 
either bosonic or fermionic. In particular, the two limits 
of the sine-Gordon model corresponding to K$g "C 1, 
where the system maps to a set of independent harmonic 
oscillators, and Ksg — 1) where the system is equivalent 
to non-interacting fermions, are examples of critical sys- 
tems with bosonic and fermionic statistics respectively. 
The other two examples in this paper (harmonic chain 
and transverse field Ising model) obviously also fall into 
the category of a system with well defined statistics of 
quasiparticles (bosonic and fermionic respectively). 

For bosonic excitations it is straightforward to show 
that the Gaussian ansatz (14311 still holds at finite tem- 



peratures for the Wigner function [3l|, [68J. The width 
of the Wigner function satisfies a similar equation as the 
width of the wave function (gH), QIOIP - The only differ- 
ence with the initial ground state is that this width gets 
"dressed" as at coth(e^./2T) at initial time, where <r° is 
the ground state width defined earlier and e° is the initial 
energy of the particle with momentum fc. Using this fact 
it is straightforward to show [3l| that: 



T cir 



coth I ^ 



Since we are interested in the excitations created by 
the dynamical process, we need to subtract from this 
quantity the initial number of quasiparticles excitations, 
which were present due to initial thermal fluctuations: 



T cff 



^Ih- coth f — 



coth I £ 



r 



n c ° x (fc)coth(^ 



, (120) 



where n|! x (fe) is the number of created quasiparticles with 
momentum k at zero temperature. Note that this result 
does not depend on the details of the quench process, i.e. 
whether it is fast or slow. At low temperatures T <C e° 
the expression (|120p obviously reduces to the zero tem- 
perature result. However, at high temperatures T ^> (° k 
we have the bosonic enhancement of the transitions 



(121) 



When we sum n^ x (k) over all possible momenta k starting 
from a gapless system, this extra T/e" factor makes the 
result more infrared divergent. In the regime of validity of 
the adiabatic perturbation theory it changes the scaling 
for example for the density of excitations to 



T\6\ 



(d-z)v/{zv+l) 



(122) 



The crossover, where the quadratic scaling is restored, is 
determined by the equation (d — z)vj{zv + 1) = 2. We 
note, however, that in small dimensions the response of 
the system can become non-adiabatic [3l[ and the scaling 
above can break down. For the sine-Gordon model in the 
bosonic limit using Eq. (|104p it is indeed straightforward 
to see that 



T(cS) 1/3 L, 



(123) 



where L is the system size. If we consider massive bosonic 
theories in higher dimensions we would find that the scal- 
ing (|122[) would be restored above two dimensions. Even 
though in one dimension the adiabatic perturbation the- 
ory fails to predict the correct scaling, it unambiguously 
shows that the finite temperatures make the response of 
the system less adiabatic. For sudden quenches the effect 
of initial temperature on harmonic system was recently 
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considered in Ref. [69( and the authors obtained results 
consistent with the statements above. 

The scenario becomes quite opposite in the fermionic 
case, Ksg — 1) where we are dealing with a sum of in- 
dependent two-level systems. At finite temperature each 
level is occupied according to the Fermi distribution 



Ik = ex P 



T 



1 



The probability of excitation thus gets corrected as: 
««(*) = <(k)(fk - ft) = »«(*) tanh Q|) . (124) 

The fact that we got tanh factor for fermions (and coth 
for bosons) is hardly surprising. Similar factors appear in 
conventional fluctuation-dissipation relations [701 ] . How- 
ever, this mere change of coth to tanh factor has an im- 
portant implication. As in the case of bosons, in the 
small temperature limit this additional factor reduces to 
unity and the zero temperature result is recovered. At 
high temperatures T 3> e° we find 



(125) 



Therefore the number of created quasiparticles is much 
smaller than in the zero temperature case. This fact re- 
flects fermionic anti-bunching. In other words the preex- 
isting thermal quasiparticles are blocking the transition 
to the already occupied excited states and the dynamical 
process becomes more adiabatic. Thus for the density of 
excitations we expect now the scaling 

n cx cx \5\id+z>/{z»+l)j T , ( 12 6) 

This analysis is not sensitive to the details of the pro- 
cess, hence the situation remains the same for sudden 
quenches. Thus instead of the scaling n cx ~ |A*| dly we 
will get n ox ~ \\*\{dTz)u T ±i wliei . e the upper « + » or «_» 

sign corresponds to bosons and the lower sign does to 
fermions. The crossover to the quadratic linear response 
scaling happens when the corresponding exponent be- 
comes two. As in the case of slow quenches one should 
be careful with the validity of the perturbative scaling in 
low dimensions for bosonic excitations, where the system 
size can affect the scaling and change the exponent. 

The main conclusion of this section is that at finite 
initial temperatures the statistic of quasiparticles qual- 
itatively changes the response of the quantum critical 
system to quenches (independently if they are slow or 
fast). At zero temperature the quasiparticle statistics 
does not seem to play an important role, since it does 
not enter the general scaling exponents. This suggests 
that the non-adiabatic response will be extremely inter- 
esting in the systems with fractional statistics of excita- 
tions. This sensitivity and potential universality of the 
response of the system to sudden or slow perturbations at 
finite temperatures might allow one to use non-adiabatic 
transitions as an experimental probe of the quasiparticle 
statistics. 



VII. CONCLUSIONS 

In this work we focused on the analysis of the response 
of a translationally invariant system, initially prepared 
in the ground state, to linear quenches, where an ex- 
ternal parameter globally coupled to the whole system 
linearly changes in time. Using the adiabatic pertur- 
bation theory we showed how to obtain the scaling of 
various quantities like the density of quasiparticles n ex , 
heating (excess non-adiabatic energy) Q and entropy Sd 
with the quench rate S, for small S. We started from a 
simple two-level system, where we showed that the tran- 
sition probability scales quadratically with 6, at 8 — * 0, 
if the external parameter changes in the finite range. We 
then showed that this quadratic scaling can be violated 
in various low-dimensional systems, especially quenched 
through singularities like quantum critical points. This 
violation comes because excitations of low energy lev- 
els, for which dynamics is diabatic, dominate the scaling 
of n cx and other quantities, which then acquire univer- 
sal non-analytic dependence on the rate 6. For example, 
we argued that in generic low-dimensional gapless sys- 
tems n e x oc \5\ d / z as long as d/z < 2. In the opposite 
limit d/z > 2 we expect that the high energy quasiparti- 
cles dominate n cx and the quadratic scaling is restored: 
n cx oc S 2 . Similar story is true for heat and entropy. We 
note that the quadratic scaling of the transition probabil- 
ity to highly excited states is specific to linear quenches, 
where the time derivative of the tuning parameter has 
a discontinuity at initial and final times. E.g. we ana- 
lyzed the situation where A = 6 for t > ti and A = 
for t < ti. If the tuning parameter is turned on and 
off smoothly in time then this quadratic scaling is no 
longer valid. One can argue, however, that at sufficiently 
high energies there are always dephasing mechanisms for 
quasiparticles which reset their phase and effectively re- 
set the value of ij . This resetting of the phase will likely 
restore the quadratic scaling. However, this issue needs 
to be investigated separately and is beyond the scope of 
this work. 

Using the same adiabatic perturbation theory and the 
general scaling arguments we showed how the universal 
behavior of n ox , Sd, and Q emerges for quenching through 
quantum critical points. In particular, we showed that 
the density of generated quasiparticles (and entropy) 
scale as n cx ,Sd ~ \5\ dl/ '( zu+1 ' , where v is the critical 
exponent characterizing the divergence of the correlation 
length near the phase transition in agreement with earlier 
works [27], . Based on the scaling analysis we showed 
that for linear quenches this scaling is valid as long as 
dv/(zv + 1) < 2, otherwise the quadratic scaling is re- 
stored. We also discussed the connection between adia- 
batic and sudden quenches near a quantum critical point. 
We argued that if dv < 2 the two are qualitatively similar 
provided that one correctly associates the quench ampli- 
tude A* and the quench rate 5: A* ~ |£|V(*"+i) . Us- 
ing the adiabatic perturbation theory adopted to sudden 
quenches of small amplitude we showed that the den- 
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sity of excitations scales with the quench amplitude as 
n GX oc |A*| d,y . When dv > 2 the quadratic scaling, which 
follows from the conventional perturbation theory (linear 
response), is restored and n cx ~ (A*) 2 . As in the case 
of slow quenches the quadratic scaling comes from the 
dominant contribution of the high energy quasiparticles 
to n ox . 

We also discussed the situation where a gapless system 
is initially prepared at finite temperature. In this case we 
argued that the statistics of the quasiparticles becomes 
crucial. In particular, for bosonic quasiparticles the fi- 
nite temperature enhances the effects of non-adiabaticity 
due to bunching effect, while conversely for fermions the 
response becomes more adiabatic than at zero tempera- 
ture due to the Pauli blocking. In the regime of validity of 
the adiabatic perturbation theory we argued that in the 
scaling laws presented above (both for sudden and slow 
quenches) one should change d — ► d—z in the bosonic case 
and d — > d + z in the fermionic case. We note, however, 
that in the bosonic case the adiabatic perturbation the- 
ory can breakdown due to overpopulation of low energy 
bosonic modes and the system can enter a non-adiabatic 
regime [3l|. We illustrated our general statements with 
explicit results for various solvable models. 
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APPENDIX A: HALF LANDAU-ZENER 
PROBLEM 

Here we briefly describe the derivation of the transition 
probability in the LZ problem, when the system starts 
at the symmetric point with the smallest gap. Let us 
consider the Hamiltonian: 



where tan0(i) = g/\(t), corresponding to the eigenergies 



H 



Xt) 9 
9 -A(i) 



(Al) 



where g is a constant and A(i) = St. This Hamiltonian 
has normalized eigenvectors given by Eq. 



sin(0/2) \ 
-cos(0/2) J 



1+) 



cos(0/2) \ 
sin(0/2) J 



(A2) 



(A3) 



The vectors representing the ground state |— ) and the 
excited state |+) significantly change as the time t is 
varied, as it is sketched in Figl3] 



i 
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energy 

^^^^^^ 











l-H(S) 
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FIG. 3: Landau-Zener problem: the eigenvectors of the sys- 
tem are exchanged during the time evolution. 



Unlike in Sec. Ill Al for the exact solution it is convenient 
to work in the fixed basis: 



hK*)> = &(*) (J) +&(*)(?) 



Then the equations of motion for the coefficients <fii(t) 
and 02 (t) become 



i<j>i = Stfa + g4> 2 
i4>2 = 94>\ - St<p2- 



(A4) 
(A5) 



This system of equations has the following generic solu- 
tions for 0i and <p2- 



4>i (t) = cje" 
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where M(a, b, z) is the confluent hypergeometric function 
otherwise called \F\. The coefflcients c\, 0%, C3, C4 are de- 
termined by the initial conditions on the wave function 
at the initial time t-,„ 



\Htin)) =M*ir 



faith 



plus two auxiliary conditions e.g. 

■ 1 I 9 , 

ifa U=o = 



lfa\t=0 



fj 



•>2t=Q 



n\t=o- 



(A8) 
(A9) 



The conventional LZ problem [46|, |47j considers the case 
of starting in the ground state at t = — 00 and asks what 
is the excitation probability after evolving the system to 
a final time t — +00, which is 



exp 



S 



(A10) 



This result can be recovered imposing the initial condi- 
tions: 



-00 |^(-oo)) = |-> 



fa (-00) = 1 
^ 2 (-oo) = 



and looking at the asymptotic behavior of |(^i(oo)| 2 . 

Another situation that can be straightforwardly con- 
sidered corresponds to starting from t = (instead of 



t = —00) and ending at t — +00, essentially this is 
half of the usual LZ-problem. As we have explained in 
Sec. IIVB~2l this situation naturally arises in e.g. describ- 
ing the problem of loading hard-core bosons into an opti- 
cal lattice. Here we need to impose the initial conditions 



Un = |V(0)> = |-} = 



V2 \ 1 



H0) = rh • 



The probability of excitation = lim t _ 
turns out to be (see also Refs. [541 l65l]) 
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I0iMI 2 
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where we introduced t = 2g 2 /S. This function has the 
following asymptotic behavior: for r — > 



and for r — > co 



1 

IT 2 " 



I65 



4 ' 



As expected the latter slow asymptotic agrees with the 
result of the adiabatic perturbation theory 



[1] I. Bloch, J. Dalibard, and W. Zwerger, Reviews of Mod- 
ern Physics 80, 885 (2008). 

[2] M. Greiner, O. Mandel, T. W. Hansen, and I. Bloch, 
Nature 419, 51 (2002). 

[3] C. D. Fertig et al., Phys. Rev. Lett. 94, 120403 (2005). 

[4] A. K. Tuchman, C. Orzel, A. Polkovnikov, and M. Kase- 
vich, Phys. Rev. A 74, 051601 (2006). 

[5] T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 
900 (2006). 

[6] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, 
and D. M. Stamper-Kurn, Nature 443, 312 (2006). 

[7] S. Hofferberth, I. Lesanovsky, B. Fischer, T. Schumm, 
and J. Schmiedmayer, Nature 449, 324 (2007). 

[8] S. Trotzky et al., Science 319, 295 (2008). 

[9] M. Rigol, V. Dunjko, and M. Olshanii, Nature 452, 854 
(2008). 

[10] P. Reimann, Phys. Rev. Lett. 101, 190403 (2008). 
[11] A. Silva, Phys. Rev. Lett. 101, 236803 (2008). 
[12] M. Rigol, arXiv:0904.3746 (unpublished). 
[13] E. Altman and A. Auerbach, Phys. Rev. Lett. 89, 250404 
(2002). 

[14] A. Polkovnikov, S. Sachdev, and S. M. Girvin, Phys. Rev. 

A 66, 053607 (2002). 
[15] R. A. Barankov, L. S. Levitov, and B. Z. Spivak, Phys. 

Rev. Lett. 93, 160401 (2004). 
[16] P. Calabrese and J. Cardy, Phys. Rev. Lett. 96, 136801 



(2006). 

[17] P. Calabrese and J. Cardy, J. Stat. Mech: Th. and Exp. 

P06008 (2007). 
[18] K. Sengupta, S. Powell, and S. Sachdev, Phys. Rev. A 

69, 053616 (2004). 
[19] C. Kollath, A. M. Lauchli, and E. Altman, Phys. Rev. 

Lett. 98, 180601 (2007). 
[20] E. A. Yuzbashyan, B. L. Altshuler, V. B. Kuznetsov, and 

V. Z. Enolskii,' Phys. Rev. B 72, 220503(R) (2005). 
[21] G. Roux, Phys. Rev. A 79, 021608 (2009). 
[22] V. Gritsev, E. Demler, E. Lukin, and A. Polkovnikov, 

Phys. Rev. Lett. 99, 200404 (2007). 
[23] S. R. Manmana, S. Wessel, R. M. Noack, and A. Mura- 

matsu, Phys. Rev. Lett. 98, 210405 (2007). 
[24] L. Mathey and A. Polkovnikov, arXiv:0904.2881 (unpub- 
lished) . 

[25] A. Iucci and M. A. Cazalilla, arXiv:0903.1205 (unpub- 
lished) . 

[26] P. Barmettler, M. Punk, V. Gritsev, E. Demler, and 
E. Altman, Phys. Rev. Lett. 102, 130603 (2009). 

[27] A. Polkovnikov, Phys. Rev. B 72, R161201 (2005). 

[28] W. H. Zurek, U. Dorner, and P. Zoller, Phys. Rev. Lett. 
95, 105701 (2005). 

[29] J. Dziarmaga, Phys. Rev. Lett. 95, 245701 (2005). 

[30] R. W. Cherng and L. S. Levitov, Phys. Rev. A 73, 043614 
(2006). 



20 



[31] A. Polkovnikov and V. Gritsev, Nature Phys. 4, 477 
(2008). 

[32] A. Altland and V. Gurarie, Phys. Rev. Lett. 100, 063602 
(2008). 

[33] R. Barankov and A. Polkovnikov, Phys. Rev. Lett. 101, 
076801 (2008). 

[34] R. Bistritzer and E. Altman, PNAS 104, 9955 (2007). 
[35] C. De Grandi, R. A. Barankov, and A. Polkovnikov, 

Phys. Rev. Lett. 101, 230402 (2008). 
[36] K. Sengupta, D. Sen, and S. Mondal, Phys. Rev. Lett. 

100, 077204 (2008). 

[37] D. Sen, K. Sengupta, and S. Mondal, Phys. Rev. Lett. 

101, 016806 (2008). 

[38] A. Altland, V. Gurarie, T. Kriecherbauer, and 
A. Polkovnikov, Phys. Rev. A 79, 042703 (2009). 

[39] U. Divakaran, V. Mukherjee, A. Dutta, and D. Sen, J. 
Stat. Mech. p. P02007 (2009). 

[40] D. Chowdhury, U. Divakaran, and A. Dutta, 
arXiv:0906.1161 (unpublished). 

[41] K. Sengupta and D. Sen, arXiv:0904.1059 (unpublished). 

[42] A. Itin and P. Torma, preprint: arXiv:0901.4778 (unpub- 
lished). 

[43] D. Rossini, A. Silva, G. Mussardo, and G. E. Santoro, 

Phys. Rev. Lett. 102, 127204 (2009). 
[44] F. Pollmann, S. Mukerjee, A. G. Green, and J. E. Moore, 

arXiv:0907.3206 (unpublished). 
[45] C. De Grandi, V. Gritsev, A. Polkovnikov, 

arXiv:0909.5181, arXiv:0910.0876 (unpublished). 
[46] L. Landau, Phys. Z. Sowj. 2, 46 (1932). 
[47] C. Zener, Proc. R. Soc. 137, 696 (1932). 
[48] T. W. B. Kibble, J. Phys. A 9, 1387 (1976). 
[49] W. H. Zurek, Phys. Rep. 276, 177 (1996). 
[50] G. Rigolin, G. Ortiz, and V. H. Ponce, Phys. Rev. A 78, 

052508 (2008). 

[51] L. Landau and E. Lifshitz, Quantum Mechanics: Non- 
Relativistic Theory, Vol. 3 (Butterworth-Heinemann, 
1981). 

[52] R. Shankar, Principles of Quantum Mechanics (Springer, 

New York, 1994). 
[53] N. V. Vitanov and B. M. Garraway, Phys. Rev. A 53, 



4288 (1996). 

[54] N. V. Vitanov, Phys. Rev. A 59, 988 (1999). 

[55] B. L. Gelmont, V. I. Ivanovomskii, and I. M. 

Tsidilkovskii, Uspekhi Fizicheskikh Nauk 120, 337 

(1976). 

[56] A. H. C. Neto, F. Guinea, N. M. R. Peres, and K. S. N. 

andA. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 
[57] T. Giamarchi, Quantum Physics in One Dimension 

(Clarendon Press, Oxford, 2004). 
[58] S. Sachdev, Quantum Phase Transitions (Cambridge 

University Press, Cambridge, 1999). 
[59] M. Girardeau, J. Math .Phys. 1, 516 (1960); C. N. Yang 

and Y. P. Yang, J. Math .Phys. 10, 1115 (1969); 

L. Tonks, Phys. Rev. 50, 955 (1936). 
[60] B. Paredes, A. Widera, V. Murg, O. Mandel, S. Foiling, 

I. Cirac, G. V. Shlyapnikov, T. W. Hansch, and I. Bloch, 

Nature 429, 277 (2004). 
[61] A. Polkovnikov, Phys. Rev. Lett. 101, 220402 (2008). 
[62] R. Barankov and A. Polkovnikov, arXiv:0806.2862 (un- 
published) . 

[63] V. Gritsev, A. Polkovnikov, and E. Demler, Phys. Rev. 

B 75, 174511 (2007). 
[64] H. Biichler, G. Blatter, and W. Zwerger, Phys. Rev. Lett. 

90, 130401 (2003). 
[65] B. Damski and W. H. Zurek, Phys. Rev. A 73, 063405 

(2006). 

[66] E. Altman and A. Vishwanath, Phys. Rev. Lett. 95, 
110404 (2005). 

[67] S.-J. Gu and H.-Q. Lin, Europhys. Lett. 87, 10003 
(2009). 

[68] A. Polkovnikov, arXiv:0905.3384 (unpublished). 
[69] S. Sotiriadis, P. Calabrese, and J. Cardy, arXiv:0903.0895 
(unpublished) . 

[70] A. A. Abrikosov, L. P. Gor'kov, and I. E. Dzyaloshin- 
ski, Methods of quantum field theory in statistical physics 
(edited by R. A. Silverman, Courier Dover Publications, 
1975). 

[71] We note that there is a small error in Ref. [27] . which 
gives a different expression for d c . 



